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■ We present an effective action approach for the problem of Coulomb blocking of tunneling. The 
, method is applied to the "strong coupling" problem arising near zero bias, where perturbation theory 

diverges. By a semiclassical argument, we obtain electrodynamics in imaginary time, and express the 
I anomaly through exact conductivity of the system a{uj, q) and exact interaction. The calculation is 
\^ ' tested by comparing with the known perturbation theory result for the diffusive anomaly. Also, we 
use the method to study the anomaly enhancement due to external magnetic field, and the effect of 
, screening by electrodes. 

a 

^ . I. INTRODUCTION 

o 
o 

Suppression or enhancement of the tunneling conductivity near zero bias is known to be a 
signature of interaction in the system. It is called "zero-bias anomaly" , and it has been studied 



;_j ■ in metals and semiconductors since the early sixties Initially, the origin of the anomaly was 



attributed to the Kondo effect due to magnetic impurities. However, later it was understood 
that a much more common mechanism is Coulomb blocking of the tunneling. A perturbation 
theory of this effect was developed by Altshuler, Aronov and Lee 0. The theory deals with the 
diffusive limit, and shows that the blocking increases at small bias which leads to a singularity 
in the tunneling conductivity. The theory has been thoroughly tested experimentally and found 
to be extremely accurate when the anomaly is a weak feature on the top of a large constant 
conductivity. 

In the last years the interest shifted to the systems with strong Coulomb effects, such as 
disordered metals and semiconductors near metal- insulator transition |^. It is found that 
the Coulomb anomaly is sharply enhanced near the transition, providing a test of Coulomb 
correlations. Another important discovery is the observation by Ashoori et al. of the Coulomb 
blocking of tunneling in a two-dimensional metal in magnetic field In this experiment it is 
found that at certain magnetic field the zero-bias anomaly abruptly increases and transforms 
to a "soft Coulomb gap" . It has been pointed out [§] that this transition is induced by disorder. 
More recently, the gap was studied in the systems with higher mobility where current is 

almost entirely blocked below certain threshold bias. These findings caused a lot of theoretical 
work, concerned with the behavior of current near zero bias |Tl|], and with determining the gap 



width |jl2|]. The anomaly has been shown to be particularly interesting in the z/ = 1/2 Quantum 



Hall system fpJll as a sensitive probe of the Quantum Hall state. 



Our goal in this contribution is to describe an effective action theory that treats the 
anomaly as cooperative tunneling. In the calculation, the tunneling rate is related to the action 
for charge spreading and is expressed through the actual conductivity cr{uj, q) of the system. 
The treatment is non-perturbative, and remains accurate in the strong coupling regime. For 
example, at the metal-insulator transition, conductiviy has scale invariant form (j{uj) ~ cj", 
and we are able to predict the form of the tunneling I — V curve. This relation may be used 
to determine the critical exponent a. We discuss it below, and also apply the method to 
diffusive anomaly in a two dimensional metal, with and without magnetic field. Our results are 



consistent with the perturbation theory results [Q,|Tl 



II. QUALITATIVE DISCUSSION 

Tunneling of an electron into a metal involves two steps: traversing the barrier, followed 
by spreading within the metal. Typically, the traversal time is much shorter than the relax- 
ation time in the electron liquid. Accordingly, it is legitimate to separate the tunneling into a 
single-electron and many-electron parts, and to treat them separately. The first contribution is 
described by the transmission coefficient of the barrier, at small bias being just a constant. We 
are interested here in the multi-electron effect that involves motion of a large number of elec- 
trons in order to accomodate the new electron. At low bias this collective effect may completely 
control the tunneling rate. 

Let us illustrate the effect of charge relaxation on the tunneling rate by using the example of 
a two-dimensional conducting plane [0 . Charge relaxation in a two-dimensional conductor is a 
classic electrodynamics problem studied by Maxwell who gave a solution in terms of a moving 
image charge |T^. In this problem, a point charge e is injected into a conducting sheet with 
conductivity cr, and one is interested in the time-dependence of the density and potential of 
spreading charge. The solution, as Maxwell formulates it, is that the potential within the sheet 
is given by that of a point charge e moving along the normal to the plane with the velocity 
V = 2na. The size of the charge cloud grows as r(t) ~ vt. Let us consider the Coulomb part of 
the action for the charge: 

^ ' J to r(t') 2n(x \toJ 

In the semiclassical picture the action (|1|) must be added to the under-barrier action. The 
divergence of (|l|) at t ^ oo indicates that for a charge spreading that takes a long time the 
spreading action dominates the tunneling rate. We will see that the time of spreading diverges 
at small bias, t^: = e/aV. From that, near zero bias the tunneling acquires a suppression factor 
of exp (^—^S{t:t)\ The estimate (|I|) showing that the action diverges at small bias means that 
the semiclassical treatment is meaningful even for a well conducting metal. However, in the 
diffusive limit the estimate (|l|) does not agree with the perturbation theory. We shall see that 
the reason is that the main part of the action is rather Ohmic than Coulomb, and that after 
writing the action properly the semiclassical method completely recovers the pertubation theory 
result. 

The relevance of the semiclassical picture in this problem can be justified by a more general 
argument, not involving specific features of the charge relaxation in two dimensions. Let us 
consider a situation when at small bias one electron crosses the barrier. Since the barrier 
traversal time is much shorter than the relaxation time in the metal, while the electron traverses 
the barrier other electrons practically do not move. Thus instantly a large electrostatic potential 
is formed, both due to the tunneling electron itself, and due to the screening hole left behind. 



(1) 



The jump in electrostatic energy by an amount much bigger than the bias eV means that right 
after the one electron transfer we find the system in a classically forbidden state "under" the 
Coulomb barrier. In order to accomplish tunneling, the charge yet has to spread over a large 
area, so that the potential of the charge fluctuation is reduced below eV. If the conductivity is 
finite, the spreading over large distance takes long time, and thus the action of this cooperative 
under-barrier motion is much bigger than h. 



III. THE ACTION FOR LONG- WAVELENGTH MODES 

For the electrodynamics problem the action can be written in terms of charge and current 
densities p{r,t) and t). Full action would also contain electromagnetic potentials, but in 
the quasistationary limit, c — oo, which we always assume below, the potentials are "slaved" 
to charges, and thus can be integrated out. As it was argued above, the contribution to the 
action of the spreading charge is mainly coming from long times when the charge deviation 
from equilibrium is small. Therefore, we can expand the action in powers of p{r,t) and j{r,t), 
and keep only quadratic terms. The action should reproduce the classical electrodynamics 
equations: the Ohm's law and charge continuity. 

Of course this requirement is entirely sufficient to determine the form of the action. However, 
it is more convenient to argue in the following way. We are going to use the action to study 
the dynamics in imaginary time. Therefore, the action is precisely the one that appears in the 
quantum partition function. The latter action expanded up to quadratic terms in charge and 
current density must yield correct Nyquist spectrum of current fluctuations in equilibrium: 



((g",gg-c.,-g)) = C^a/3kl + aa' D pp: q^' Qp' . (2) 

Here 

g=j + DVp (3) 

is external current and Dap is the tensor of diffusion constants related to the conductivity 
tensor by the Einstein's formula: a = e^uD, where u = dn/dp is compressibility. Generally, 
both a and D are functions of the frequency and momentum. For simplicity, we assume that 
the temperature is zero and discuss only a two dimensional metal with spatially isotropic and 
homogeneous conductivity: axx = CTyy, a^y = —(Tyx- 

The requirement that the action produces correct current fluctuations is essentially equiv- 
alent to the fluctuation-dissipation theorem. Thus, the form of the action is fixed by response 
functions of the system. In imaginary time we get 



S = ^ I I d'^Xid'^X2 



gl ^xi-X2g2 + 



(4) 



where xi^2 = (^1,2, ^"1,2), and 612 = S{ti — ^2)- The kernel Kr,t is related to the current-current 
correlator, 

iK-!,)ap = {{SI,S-^^,-,)) (5) 

given by (0), where a and D are functions of the Matsubara frequency related with the real 
frequency functions by the usual analytic continuation. We take Coulomb interaction in the 
second term of the action (^) as non-retarded because we are going to study systems with 
relativelv low conductivitv, and thus slow charge relaxation. 



One can also justify the form @ of the action by looking at various particular examples, like 
the ideal dissipationless liquid of charge, or the phenomenological dissipative Caldeira-Leggett 
system with spatially distributed coupling to a thermal bath. In the first case the Ohmic part 
of the action is replaced by the kinetic energy term 



'kin 



j2 

drdt 

2ep 



In the second case one has 



ScL = I drduj 



ATie'^ri\uj\ ' 

where rj is the Caldeira-Leggett viscosity. For Skin the electric impedance of the system is 
imaginary, while for Sql it is real. From these two "limiting" cases one can conjecture the 
more general form (0). 



IV. INSTANTON IN IMAGINARY TIME 

To evaluate the tunnelling rate, we use the instanton method and look for a least action 



"bounce" path in imaginary time ||T^,[T6|. Among the bounce paths symmetric in time, p(r, t) = 
p{r, —t), j(r, t) = — j(r, —t), we shall find the least action path which will give a semiclassical 
estimate of the tunneling rate exponent. 

From the variational principle one can derive equation of motion. We note that the action 
@ contains the charge and current densities as miie^eniient variables, as Eq.@ was derived by 
matching with the equilibrium fluctuations in the grand canonical ensemble where charge is not 
conserved. Therefore, we have to supply the action (|^) with the charge continuity constraint: 
p + V ■ j = J7'(r, t), where 

J{r,t) = e6{r){6{t + T)-6{t-T)) . (6) 

This form of the charge source J7'(r, t) describes electron injected the system at r = 0, t = — r, 
and taken back at t = r, at the same point. In principle, one could consider a more general 
source term J' of the form J' = e{6{t + r)5(r — ri) — 6{t — r)5(r — r2)), which corresponds to 
the process of one electron entering and then leaving the liquid at different points. However, in 
real situation the tunnelling occurs preferentially at "hot spots" , or defects, where the tunnelling 
barrier is low or narrow. Thus we assume ri = r2, which also accounts for the quantum point 
contact tunneling experiment. 

The charge continuity requirement is incorporated in the action by adding a Lagrange 
mutiplier term: 

Stotai=S{p,i)+(P{r,t){p + Vj-J{r,t)) , (7) 

where 0(r, t) is an independent variable of the problem. For the least action path, the variation 
of Stotai relative to infinitesimal change 6(f), 6p, and (5j vanishes. (Note that due to the charge 
continuity 6p + V ■ 6j = 0.) After eliminating we get the standard equations of classical 
electrodynamics: 

(i) p + V-j = J{r,t) ■ 

(ii) j + DVp = a(cu,g)E; (8) 

(iii) E(r, t) = -Vr [ dr'p{r', t)U{\r -r'\) . 



The equations describe the system trajectory in imaginary time, i.e., the spreading of charge 
under the Coulomb barrier due to selfinteraction. 

Then one must solve Eq.(|) for p and j, and to compute the action (^). For a spatially 
homogeneous system, by using Fourier transform, we get 

p[uj,q) - 



\uj\ + Dq^ + cr^^q'^Ug 
}{^,<l) = -iK-\uj,q)qUgp{uj,q) , 

where Uq is the Coulomb potential formfactor. Substituted in Eq.(^) this yields the action 

S(r)-W (9) 
' 2i-Juj\ + Dq^ \uj\ + Dq^ + a^Ug ^ ' 

which depends on the accomodation time r through Fourier component of the charge source: 
J^{uj) = 2iesmujT. 

Finally, to obtain total action of the system we subtract from the action So{t) of spreading 
charge the term 2eVT that accounts for the work done by voltage source: S{t) = iSo(r) — 2eVT. 
Thus the energy conservation at transferring one electron across the barrier is assured. Then 
one has to optimize S{t) in r. Optimal satisfies the relation 

?^^ = 2eV (10) 

Having found from Eq.(|TD|) one gets the tunneling rate that coinsides with tunneling con- 
ductivity up to a constant factor: 



G{V) = Goexp 



-Uso{r.{V))-2eVT4V)) 
n 



The optimal r* can be interpreted as the charge accomodation time. 

Let us point out here that the accuracy of the term 2eVT is determined by the assumption 
that ^ Tf, the time it takes one electron to traverse the barrier. This assumption is valid 
whenever there is an anomaly: if ^ tj, then S{t^) ~ h, and thus there is no tunneling 
suppression. 

To summarize, the equations (^, ([T0|), and (pH]), taken together, define tunneling conduc- 
tivity. Within this quite general framework, one can study the anomaly in different systems. 
Let us emphasize that after having calculated t^{V) and 5(r*) it is essential to check the self- 
consistency of the assumption that r^, ^ tj. For example, this assumption will not be fulfilled 
in a clean metal, i.e., in Fermi liquid without disorder {D > 1). The reason is that the con- 
ductivity of an ideal conductor is <j{uj) = ine^ /moj. In this case Eq.(|D gives Soij) h at 
any t ^ Tf, and henceforth t^, c:^ Tf. This indicates absence of the anomaly in a clean metal, 
the result familiar from the Fermi liquid picture. On the contrary, for a one- dimensional metal 
iSo(t) ~ In t/tj, which leads to the power-law anomaly known from the Luttinger liquid theory. 



V. COMPARISON WITH DIFFUSIVE ANOMALY 

For a two-dimensional metal with elastic scattering time tq and non-screened Coulomb 
interaction we set Ug = 27i/\q\ and a^^ = a, constant at |a;|,Wi7|g| < I/tq. Then Eq.(^ gives 

•^M-sitaf;)'"!--"-^) . (12) 



(Here v is compressibility.) From Eq.(p!0|), 



= ^ X^iUovelV) . (13) 

The theory is selfconsistent in the hydrodynamic hmit, r* > tq, i.e., at eV < e^/aTQ. Then the 
least action is 

/ e \ 1 / ero(T(i^e^''^^ 



5(1^) = — - In —— In n^^l^i:,^ (14) 

It is interesting to compare this result with the identical double-log dependence derived by Alt- 
shuler, Aronov, and Lee in a different context 0. They calculated perturbatively the correction 
to the tunneling density of states 5v{t) with the assumption that it is small, \5v\ <^ z/q, which 
is the case only for a weak disorder. It was found that 5z/(e) = —Ti^^u^SiV = e/e), where S(y) 
is given by ([T^). The main difference is that our double-log has to be exponentiated to get 
the tunnelling density of states, while in the double-log itself appears as a correction to the 
density of states. In the range of the perturbation theory validity the two results agree. From 
that point of view, our calculation provides description of the diffusive anomaly at low bias, 
where the perturbation theory diverges. 



VI. SCREENING BY ELECTRODES 



In a real experiment the charge tunnels between two electrodes, and there are separate 
contributions to the action due to the relaxation of the electron and hole charges on both sides 
of the barrier. If the electrodes are close, the charges partially screen the field of each other, 
which makes their spreading correlated. In this case the least action is smaller than the sum 
of independent contributions of the electrodes, and thus the anomaly is weakened. For a two 
dimensional system, quantitatively, the effect will be that the log-divergence of the integral over 
q in Eq.(P) will be cut at g ~ a~^, where a is the distance between the electrodes. As a result, 
the l^— dependence of the second log in Eq.(p!^ saturates at eV Vq = ha / a. 

This "excitonic" correlation effect can be treated straightforwardly by writing the action (^ 
for each electrode separately, together with the term describing interaction across the barrier. 
Let us consider an example of two parallel planes with different conductivities and diffusion 
constants, cxi, (T2, Di, and D2. It is straightforward to generalize the action and to find the 
instanton. The least action is 



Mr) 



u),q 



cj| +Dg2 + Sg2u) 





■ 1 " 




.-1. 



(15) 



where we use matrix notation: 
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iqr 



(16) 



The rows and columns of the matrices (|T6D correspond to the two planes. 

Let us assume e^i/i^2'2 ^ 1, which is the case in almost all experiments. Then, by carrying 
out matrix inversion and integration we get 



5nfr) 



«ln (—^ at r > h/eVo . 



;i7) 



Here 



47r2 



1 Airaa 1 Anaa 
m 1 111 



(18) 



where a = o-i,xxO'2,xx/ {cri^xx + cr2,xx)- If the two planes are identical, 



a 



■ In 27ie ua . 



(19) 



Thus at V < Vq the I — V curve becomes the power law / ~ The tunneling suppression 

in this case is weaker than for the non-screened interaction. 

Like other interaction effects, the anomaly is enhanced if the system dimension is lowered. 
Let us consider tunneling into a disordered wire of thickness d from a well conducting electrode, 
two- or three-dimensional, parallel to the wire at a distance a ^ d away from it. In this case 
the main contribution to the action is due to charge relaxation in the wire. We still can use 
Eq.(P), where J2 ■■■ = I p---- and the two-dimensional a^x is replaced by Aa, where A is the wire 

crossection and a is the three-dimensional conductivity. The Coulomb potential formfactor at 
small k -C a~^, using electrostatic image in the electrode, is found to be Uq = 2\na/d. By 
doing the sum in (H) over k and u we get 



So{t) = BV^, B 



(20) 



Then Eq.(|TUp gives optimal t^{V) = B /16e V , and thus the conductivity 

/ B^ \ 



G{V) ~ exp 



8heV 



(21) 



The "field+charge" diffusion constant alloA that appears in the ID plasmon dispersion relation 
uj{k) = —iaUoAk"^ is typically much bigger than D. Then the spreading of charge is effectively 
one dimensional at V < Vi = e'^Uo/Sa. For a thick wire Eg . (pl|) holds in a wide voltage range 
Vi > V > Vdiff where one can ignore localization corrections to cr. By the Thouless criterion, 
the voltage Vdiff above which the diffusive treatment is valid is estimated by comparing the 
tunneling time t^{V) with the inverse level spacing in the region over which the charge spreads, 
-1 



(22) 



A-i = uAJaUoAr^iV). This gives 



On the time scale T^,{Vdiff), the transport is diffusive if the wire length L is much shorter than 

huA'^aUo. 



^diff = \/<yUoAT^{Vdiff 



VII. TUNNELING IN A FINITE SYSTEM 

It is of interest to consider the tunneling problem for a system of small size in which the 
time of charge spreading is limited by the dimension. Most important realization is tunneling in 
small Coulomb blockade devices, like quantum dots or very small capacitors. For such systems, 
at sufficiently low bias, the dynamics of charge spreading is unimportant and the effect of charge 
relaxation can be taken into account by modelling the conducting part of the system by a single 
effective resistor. 



In this limit, the problem can be treated by the recently developed theory of the "environ- 



ment effect" on the tunneling conductivity P-rO|. In the latter the environment of a tunnel 



junction is replaced by an effective electric circuit that causes voltage fluctuations across the 
junction. Assuming that the fluctuation spectral density is equilibrium, one calculates the 
conductivity suppression factor by doing gaussian average of the tunneling current, in a way 
similar to the Debye- Waller theory of the fluctuation effect on the structure factor of crystals 
[^,10]|. Alternatively, one can calculate the tunneling transition rate that includes the shake up 



effect of the plasmon modes in the environment ^ . The results of Refs. |8|-[10| can be derived 
by our method if one takes the zero wavelength limit of the action (^) by keeping only the 
largest scale electromagentic mode for each element of the circuit. (For example, in a resistor 
one takes uniform current distribution and zero charge density, or in a capacitor one assumes 
uniform charge distribution over the capacitor plates and no current.) Then the instanton cal- 
culation will lead to the results identical to those of Refs. P|-pi][|. This zero wavelength limit 
describes tunneling at very low voltage, when during the under-barrier motion the tunneling 
charge spreads over a scale comparable to the size of the system. 

Therefore, there are several regimes corresponding to different values of the bias. Let us 
list them here for the problem of tunneling between two identical conducting sheets L x L of 
conductivity a, separated by the barrier of thickness a. 

• Effectively infinite system: e^/a > eV > e^a/L? 

G = dI/dVr.V-^ ^^e^ln(eVa) 



Finite system characterized by impedance of environment: e^a/L"^ > eV > hD/L 



2 



a = -—\nL/d; (24) 



• Ohm's law restored at eV < hD/LP". G=const 

In the expression for a given for the infinite system it is assumed that e^va ^ 1. The result 
for finite systems agrees with Refs. fl^TO|. The third case describes the limit of very low bias, 
when the system is characterized by an overall transmission coefficient, which is equivalent to 
connecting in series the barrier and the metal effective resistors. The three cases describe the 
change of the conductivity over the whole voltage range, from relatively high voltage where 
perturbation theory works, through the region of low voltage where the anomaly disappears. 



VIII. 2D ELECTRON GAS IN MAGNETIC FIELD 

It is straightforward to incorporate magnetic field in the theory by substituting a^^ = a^xiB) 
in the action (|T2|). (Note that a^y does not enter.) As magnetic field increases, the conductivity 
drops, and at certain field it reaches the quantum limit cxg = e'^/h. In this field range the 
prefactor a in Eq.(|l7|) becomes of the order of one, and the anomaly in the conductivity 
changes from weak to strong. The threshold conductivity, according to Eq. ([l^) , is 

1 

(Tc = — -— ln27re^i/a (25) 
4:71^ h 



A transition like that was observed by Ashoori et al. 0] in the tunneling current from a 3D 
metal into a 2D electron eas. In this exDeriment. the ohmic conductance was measured as 



function of temperature, which corresponds to our zero temperature non-hnear current taken 
at \^ ~ kBT/e. The 2D gas was relatively clean with the zero field mobility corresponding 
to the elastic scattering time tq ^ 4 ■ 10~^^ s. The Fermi energy calculated from the electron 
density was Ep ~ 10 mV. By using the result (^) together with the Drude-Lorentz model, 
o'xxiB) = ne^ /mTQUil at tquJc ^ 1, one finds that the anomaly hardening transition corresponds 
to the cyclotron frequency u* = (SnEp/hrQY^'^ ~ 8.0 mV. In terms of the field intensity this 
is approximately 4.6 Tesla which is quite close to the transition field reported in Ref |Q. 

It is interesting to note that in a weakly disordered metal with EpT^ ^ % the threshold 
field is small: UuJq <^ Ep- This means that the transition occurs well below the field where 
the Quantum Hall state is formed. Therefore, our estimate of uJq based on the "bare" Drude- 
Lorentz conductivity is meaningful and legitimate. On the other hand, to find the current at 
very low V one would have to use the conductivity renormalized by localization and interaction 
effects. 

Finally, let us mention a relation with the work by Halperin, He, and Platzman [llj that 



deals with the anomaly in the ly = 1/2 Quantum Hall state. In this work, the problem was 
treated by summing linked cluster terms of perturbation theory, with the density responce 
function borrowed from the Chern-Simons Fermi liquid theory . The anomaly was found to 
have the form: 

G{V) ~ exp i-Vo/V) , Vo = Att^^/^ , (26) 
where ^ <^ Vq, and n is density. It is interesting to see how this result can be derived from 



the effective action. It has been shown |[T^ that conductivity of the v = 1/2 state has strong 



spatial dispersion: 0"^ = A\k\, A = (P' /l^Tie^/wn. If this form is inserted in the action (||), one 
gets S{t) = ■7tJ2t/A, which leads to the tunneling rate (|26|) . 



IX. CONCLUSION 

The essential feature of the presented approach is that it relates the tunneling anomaly with 
exact conductivity of the system without calculating it. This would allow a comparison with 
experiment in the situations where there is no accepted model for conductivity. For example, 
the tunneling current as function of voltage can be taken from experiment, and directly used 
to find S{t) by inverse Legendre transformation. Thus obtained function S{t) can be analyzed 
to extract the conductivity frequency and wavevector dependence. 

To summarize, we argued that the theory of the Coulomb anomaly in the regime of strong 
suppression of tunneling is semiclassical. The underlying reason is that the transfer of one 
electron across the barrier is controlled by cooperative motion of many other electrons. We 
treat this motion as classical electrodynamics in imaginary time, write down the action and find 
the least action instanton trajectory. The results are compared to several known perturbation 
theories and an agreement is found. The theory is used to interpret experiments on the coulomb 
gap induced by magnetic field. 
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BLOCKAGE COULOMBIEN DE L'EFFET TUNNEL: 
DE L'ANOMALIE A ZERO VOLTAGE AU GAP COULOMBIEN 



Nous presentons une approche de la theorie effective pour le probleme du blockage coulombien 
dans I'effet tunnel. La methode est appliquee au probleme de la limite de fort couplage proche du 
voltage zero, ou la theorie de la perturbation diverge. Par un argument semiclassique, nous obtenons 
une theorie electrodynamique en temps imaginaire, et nous exprimons I'anomalie en terme de la 
conductivite exacte (t(w, q) du systeme et de I'interaction exacte. La validite de notre approche est 
teste par comparaison avec le calcul connu de I'anomalie diffusive obtenu en theorie de la perturbation. 
Nous utilisons aussi cette methode pour I'etude de I'accroissement de I'anomalie du a un champs 
magnetique externe, et a I'ecrantage des electrodes. 



